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Abstract 



Quasi one-dimensional systems are systems of particles in domains which are of infinite extent in one 
direction and of uniformly bounded size in all other directions, e.g. on a cylinder of infinite length. The 
main result proven here is that for such particle systems with Coulomb interactions and neutralizing back- 
ground, the so-called "jellium", at any temperature and at any finite-strip width there is translation symme- 
^ try breaking. This extends the previous result on Laughlin states in thin, two-dimensional strips by Jansen, 

Lieb and Seiler (2009). The structural argument which is used here bypasses the question of whether the 
translation symmetry breaking is manifest already at the level of the one particle density function. It is akin 
to that employed by Aizenman and Martin (1980) for a similar statement concerning symmetry breaking 
at all temperatures in strictly one-dimensional Coulomb systems. The extension is enabled through bounds 
OO which establish tightness of finite-volume charge fluctuations. 
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1 Introduction 

In this work we investigate symmetry breaking in classical quasi one-dimensional "jellium", that is particle 
systems with Coulomb repulsion and attractive neutralizing background (also known as "one-component 
plasma"), and in quantum systems whose states may be described by such ensembles. The particles are 
of equal charge —q, and move in domains which are of infinite extent in one direction and of uniformly 
bounded size in all other directions, e.g. cylinder or tube of infinite length and a finite, uniform, cross- 
section. The Coulomb potential is the solution to a Poisson equation with Neumann or periodic boundary 
conditions in the confined directions. It corresponds to the situation where not only the particles but also 
the electric field is confined to the tube. 

Our main result is that such systems display translational symmetry breaking in the long direction, e.g., 
the cylinder axis, which is denoted here by x. This generalizes previously known results in one and two 
dimensions. The proof is by a structural argument which is not limited to low temperatures or small tube 
cross sections. 

For the one-dimensional jellium this phenomenon was shown in 1 10, 4, 2 |. The periodicity is related to 
the possible formation of a "Wigner lattice". Roughly, the Coulomb interaction leads to the suppression of 
large scale deviations from neutrality. Each particle, as ranked by the x coordinate, fluctuates with only a 
bounded deviation from the lattice value of x corresponding to its rank. 

In two dimensions, with periodic boundary conditions in the confined direction, symmetry breaking is 
known for some special cases corresponding to even-integer values of the so-called "plasma parameter" 
Pq 2 (with /3 the inverse temperature). When f3q 2 — 2, the model is explicitly solvable and periodicity was 
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shown in 0. A proof of symmetry breaking for other even-integer values /3q 2 = 2p and sufficiently small 
values of the strip width (compared to (density) -1 ' 2 ) was given in (5), where the focus was on the Laughlin 
states in cylindrical geometry. Numerical results may be found in 1 17 1. 

The case of even-integer values of the plasma parameter is of an additional interest, and further tools 
are available for it, as it relates to Laughlin' s wave function II II which is frequently used as an approximate 
ground state for electrons in the context of the fractional quantum Hall effect. The function models the state 
of an electron gas whose filling factor (a quantity related to the electron density) is a simple fraction 1/p 
with p an integer. The function's modulus squared is proportional to the Boltzmann weight of a classical 
one-component plasma: |^| 2 oc exp(— /3U). The filling factor of an electron gas and the plasma parameter 
of a classical Coulomb system are related by fiq 2 = 2p. The solvable case fiq 2 = 2 corresponds to p = 1 
(filled lowest Landau level), which models non-interacting fermions. The integrality of p enters the proof 
in |9| in a crucial way, allowing to expand p-th powers of polynomial into monomials. 

The proof given here follows an altogether different approach, along the lines of (2] [T) . A key point is 
that due to the strong tendency of Coulomb systems to maintain bulk neutrality an interesting phenomenon 
occurs in one-dimensional and quasi one-dimensional Coulomb systems : the total charge in cylinders of 
arbitrary length is of bounded variance. Even in the infinite-volume limit the question "what is the total 
charge at points with x < xo?" has a well defined answer (denoted here by qK(x, uj)). As was pointed out 
in 1 1 1 that fact in itself implies symmetry breaking, through the long-range correlations in the values of the 
phase e' 2 ' !TK ( x,u \ or equivalently through the fractional part of the total charge below xo- This approach 
to symmetry breaking was developed in the earlier work on strictly one-dimensional Coulomb systems (2J, 
in which case qK(x, u) yields the electric field at a;. In 1 1 1 this was extended into a general criterion that 
in one dimension, "tightness" of charge fluctuations implies symmetry breaking. To make this argument 
applicable in our setting one needs to first establish the tightness estimates. These are somewhat more 
involved than in the strictly one-dimensional case. Curiously, in both [9| and here symmetry breaking is 
explained through a combination of analytical and topological arguments. However, these seem to be of 
somewhat different nature in these two works. 

The paper's outline is as follows. The model is introduced in Section[2] The main results are stated in 
Section[3] that includes the statement of Theorem 1 3 . 1 1 and its Corollary |3 . 2[ which addresses the symmetry 
breaking in the model's infinite-volume Gibbs states. The results are enabled by a pair of estimates which 
play an essential role: Theorems |3.3| and |3.4| We then turn to the mathematical framing of the particle- 
excess function K(x,ui), which it is convenient to view alternatively as a random element of a Skorokhod 
space, and as a cocycle in the sense of dynamical systems. These terms are discussed in Section [4] and 
used in Section|5]for a conditional proof of symmetry breaking, assuming the tightness bounds. The proof 
of the latter involves a different set of considerations. These are outlined in Section [6] which reviews 
the corresponding question in one dimension. Finally, the proof of the enabling bounds is spelled out in 
Sections[7] where we establish tightness of the distribution of K(x, cj) at fixed x, and in Section[8]where 
it is shown that the volume-averages of K(x, uS) tend to null. We end with a few additional remarks in 
Section[9] 

2 Coulomb interaction in quasi one-dimensional systems 
2.1 The potential function 

The quasi one-dimensional systems considered here are systems of particles which, along with the electric 
field they generate, are confined to a tubular region of the form T — R x D where D is a compact subset 
of R*, possibly with some periodic boundary conditions. The precise technical assumption is spelled out 
below. In the simplest example T is a strip whose cross-section is D = [0, W] with periodic boundary 
conditions. 
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Points on T will be denoted z — (x, y) with i£l and y £ D. For simplicity we denote the volume- 
form on D by dy and its total measure by W — J B dy. 

The Coulomb potential between two points is given by a symmetric function, V(z, z') = V(z! , z), 
which satisfies: 

- AV(z,z') = S(z-z') (1) 

Ft 2 

for — A = — #73- — An, with Ae the Laplacian on D which is taken here to be defined with either periodic 
or Neumann boundary conditions. 

Since T has locally the structure of R d , with d = 1 + k, the short distance behavior of the potential at 
interior points of D is 

{ \{d-2)C d }- 1 Aist(z,z')- {d - 2) ford / 2 
V(z,z") » r K \ . (2) 

I -(2 7 r)- 1 lndist(z,z') ford = 2 

Yet, at long distances Viz, z') behaves as a one-dimensional Coulomb potential: 

V{z,z) w -\x-x'\/{2W), (3) 

in a sense which we shall now make more explicit. 

In the example of the 2D periodic strip it is convenient to use the complex notation: z — x + iy, in 
terms of which 

V(z,z) = -(27r)" 1 log|2sinh(7r( 2 -z')/W / )l- (4) 
This function is clearly periodic in y = Im z and harmonic throughout the (periodic) strip except at z = 0, 
and it can be easily seen to satisfy {T}. At short distances V(z, z') behaves as the two-dimensional Coulumb 
potential, with logarithmic divergence, but its long distance behavior is close to that in one dimension, and 
better described by the decomposition: 

V{z,z') = -\x-x'\/(2W) + V 2 (y,y';\x-x'\) (5) 

with the correction to the linear term given by 

V 2 (y,y';\x-x'\) = -{2*)^ log |1 - e -^(\^'\+Mv-y'\)/w { _ 

which decays exponentially in \x — x'\. 

In the more general case the potential admits the eigenfunction expansion: 

V(z,z') = -\x-x'\/(2W) + ^(2^)- 1 e- | ^' l ^^(y)^(y') 

n>l 

=: -\x-x'\/{2W) + V2(y,y';\x-x'\). (7) 

in terms of the eigenfunctions of A: — Ao^niy) = E n (p n (y). 
We may now state our assumptions on D, which are: 

1 . The correspondingly periodic / Neumann Laplacian Ad is a self adjoint operator with a non-degenerate 
ground state (f(y) = W -1 ' 2 ) and compact resolvent. 

2. The Coulomb potential in T is of the form (|7} whose second term averages to zero over D: 

/ V 2 (y,y';\x-x'\)dy = 0, (8) 
Jb 

for any x, x' G R and y' 6 D. 
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3. The term V 2 admits bounds of the form 

V 2 (y,y';\x\) > -g(\x\) (9) 
and for each S > there exists c(<5) > such that 

V 2 (y,y';\x\) < cgQx\) for all \x\ > e (10) 
for g(x) a positive decreasing function on [0, oo), which satisfies the 'finite energy condition': 

x g(x) dx < oo . (11) 

The boundary conditions ensure that the electric flux lines (i.e., lines tangent to VV) do not leave the 
tube, see Fig. |2.1| 



E - field 



Figure 1 : Flux lines of the potential with Neumann boundary conditions. 

The averaging condition ([8} is equivalent to saying that the interaction of a uniform slab with a point 
charge does not depend on the charge's position within the tube's cross section. In terms of the eigenfunc- 
tion expansion this is implied by the constancy of the ground state <fio, to which all other eigenstates ip n 
are orthogonal. The decay in \x\ of the individual terms is due to the spectral gap. However, one has to 
control the sum in (17), whose terms need not be uniformly bounded (in L°°). We owe to Rupert Frank the 
comment that conditions (9| and {10) hold for compact Lipschitz domains. In M 2 , this class includes com- 
pact domains with piecewise differentiable boundary, which may exhibit discontinuities in the tangent's 
direction but no "horn singularities" of vanishing angle. For bounded domains D with smooth boundary, 
these conditions hold with g(x) of exponential decay. 

2.2 Quasi one-dimensional jellium 

Jellium in a finite segment, corresponding to 

T[L lt L 2 ] = [Li,L 2 ] x D with Li, L 2 G R , 

consists of a collection of N ~ p(L 2 — Li)W particles of charge (— q) each, with p the mean number of 
particles per volume, moving in a neutralizing background of homogeneous charge density qp. 

We shall denote a particle configuration by u = (zi, . . ■ , zat), with Zj = (xj, yj) £ 7[l 1 ,l 2 ] an d me 
labeling chosen in the increasing order of the ^-coordinates. Thus, the configuration space is La ] = 

A(N, [Li,L 2 ]) x with A(N, [Li,£ 2 ]) the simplex {x £ R N : L x < xi < ■ ■ ■ < x N < L 2 }. 

Using {7| and <j8j, one gets for the system's energy, discarding the finite self-interaction of the fixed 
background: 



U(z 1 ,...,z N ) = 2^ q V(zj,z k ) ~2^q P / V(z jt z)dz 

l<j<k<N j = l JT iLliL2 ] 



<7 2 V(z j ,2 fc ) - J2q 2 p j 

i<k<N 3=1 

?Vfe, Z + Jg 2 pE(^-H^) 2 + Const(iV,L). (12) 



l<j<k<N i=l 
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with dz — dx dy the Lebesgue measure on 7[l 1 ,l 2 ] an d L = L% — L\. 

The jellium's Gibbs equilibrium state, at the inverse temperature P = ^L, is the probability measure on 

nw t ; 

e-P u{uj) du [Ll<L2] /Z(p,N,L) (13) 

where doj[ LltL2 ] denotes the product Lebesgue measure on fir i;L La j ■ The normalizing factor Z(P, N, L) 
is the finite volume partition function. 

One may note that for D = [0, W] with the periodic boundary conditions and Pq 2 = 2p, with p an 
integer, this measure coincides with the Laughlin states, | ^ p (zi , . . . , zjv ) | 2 , which were considered in 1191 . 

Dimensionless parameters and choice of units One could, for convenience, rescale the length 
units (uniformly in all directions) so that the particle density is p = 1, and, since q affects the state only in 
the combination Pq 2 , absorb q 2 in the units of /3 thereby setting q = 1. With this choice one is still left with 
dependency on two parameters, P and W . We shall follow this choice when referring to constants which 
appear in bounds throughout the paper as C(/3, W). However, in other places we shall leave the p and q 
dependence explicit. 

An explicitly relevant length scale for us is 

A := (pW)- 1 . (14) 

This is the length of a cylindrical cell in which the mean number of particles is one. It is also the basic scale 
for the translation symmetry breaking proven here. 

It may be added that for dimensionless parameters of the model one may chose W /1//( - d_1 '/A = 
pl4/ d /( d_1 ) and the "plasma parameter", sometimes called the coupling constant, F^ = Pq 2 p ( - d ~ 2 ^ d . 



3 Statement of the main results 

3.1 Symmetry breaking 

Our main result is the following statement concerning the infinite-volume limits of the finite volume Gibbs 
equilibrium states of the quasi one-dimensional systems in the regions T[l 1 .l 2 ]- 

As is explained by the result, it is natural to take the boundaries of the finite regions as: 

Li = (-ni-9)X , L 2 = (n 2 -e)\ (15) 

with m, ri2 £ N, and 9 £ [0, 1). The corresponding Gibbs equilibrium measure for the neutral systems of 
N — ni + 712 particles in such domains will be denoted here pif^f ^ . For the purpose of convergence 
statements, these can be viewed as point processes on the common space T = 1 x D. 

Theorem 3.1. At all p > OandW > 0: 

1. For any 6 G [0, 1) and any sequence of integer pairs {(ni(j), n2{j))}, with ni(j),n2(j) — > oo, 
there is a subsequence for which /i„j'^ 9 ' converges to a limit, in the sense of convergence of proba- 
bility measures on the space of configurations in T (i.e. of point processes). 

2. Any two limiting measures which correspond to different values of 6 are mutually singular. 

3. More explicitly: there exists a measurable function on the space of configurations, such that for each 
9 e [0, 1); 

$(w) = e 27vB (16) 

almost surely with respect to each measure which is a limit of a sequence of finite-volume '9 Gibbs 
states '. 
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Furthermore, under the shifts (induced by T u : x h- > x + u) and reflections (induced by 1Z : x \— > 

-x): 

$(T u cj) = $(cj) e l2WA and $(7£w) = . (17) 

This has the following elementary consequence: 
Corollary 3.2 (Symmetry breaking). 

1. None of the limiting measures discussed above is invariant under the shifts T u , with u ^ AZ (at 
\=(W P )- r ). 

2. Except for the two cases 6 = 0,1/2 (for which 0=1 — 6 fno j i ), the limiting measures are also not 
invariant under the reflection 1Z. 

The derivation of the above statements combines three elements: "statistical mechanical" bounds on the 
distribution of the 'particle-excess function', which are presented below, dynamical systems' concepts of 
tight cocycles, and simple topological considerations. 

A fundamental role in the analysis is played by the 'particle-excess' function K(u), which expresses 
the difference (below u) of the total background charge times g _1 and the number of particles: 



K(u,ui) 







\{j < N : Xj < u}\ ifu e [L U L 2 ] 

ifue R\[Li,L 2 ] 



(18) 



with | ■ | the cardinality of the set (see Figure 2). For a point process on the line it is rather exceptional 
that such a quantity has a good infinite-volume limit (Li — > — oo, L-z — > oo). When it does, translation 
symmetry breaking follows by the general argument of [T]. 

1" r— 



Li 



K(x) + 




Figure 2: A particle configuration on [Li, L2] X D (here 
associated particle-excess function K(x). 



[— 7T, 7Tj, p 



1) along with its 



3.2 Bounds on the distribution of the particle-excess function 

The results stated in Theorem |3,l| are enabled by the following auxiliary bounds. The first bound expresses 
the tightness of the distribution of the particle-excess function. 
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Theorem 3.3 (Tightness bound). For each /3 > and W > 0, the following bound holds for all n\,n% G 
Z, 8 G [0, 1), and x G [Li, La] (with Lj defined in \\5\ ), 

^2 ({\K(x;lj)\ > 7 }) < C(f3,W)e- A ^ 2 (19) 

at some C(/3, W) < oo andA(/3, W) > 0. 

This bound can be made intuitive by noting that the Coulomb systems' energy can be presented as 
^ / K(x) 2 dx plus the short-range interaction corresponding to V2. In fact, a stronger statement is valid, 
with the exponent in \19\ replaced by —A(j3, W)^] 3 . The bound asserted in \19\ applies also to the two 
component Coulomb system. It is stated here in this weaker form in order to make the discussion of its 
implications applicable to also such systems. 

The second bound expresses the asymptotic vanishing of the translation averages of K[x, lj): 

Theorem 3.4 (The vanishing of K's volume-averages). For each f3, W, 8 > 0, there exist C — C'(/3,W,8) < 
00 and a — a(/3, W) > with which 

for any n\,ri2 £ Z, x G R, and any r > 0. 

Since K(x, w) = for x ^ [Li, L2], it suffices to verify l |20| > for all x G [Li, La] an d r G [0, L2 — x]. 
To simplify the notation, we shall sometimes write P instead of Mff^'*/ • 

Before presenting the derivation of these estimates (Sections|7]and[8|, we shall show how they imply 
symmetry breaking, i.e. give a conditional proof of Theorem |3.1| First however, let us frame the discussion 
of the function K(x; uj) within a convenient setup. 



^["1,12] 



/ K(x + u; u) du 
Jo 



4 Basic properties of the particle-excess function 
4.1 A convenient topological setup 

It is convenient to view the particle-excess function K(x; ui) as a random variable with values ranging over 
a subset of the Skorokhod space, which we denote by 5a , of functions on K which 

1. have the cadlag property (continuity from the right, and existence of limits from the left), 

2. have only integer valued discontinuities, 

3. are piecewise differentiable with a constant derivative, given by pW = 1/ A. 

It may be noted that for elements of S\ the cadlag continuity modulus wk(S), which expresses the maximal 
variation of K(x; lj) between pairs of points at distance 8 omitting a finite number of discontinuities (3), 
takes the non-fluctuating values: 

w K (8) = 8/X, for all 8 > and K G S x ■ (21) 

We shall make use of the following two observations, which are simple consequences of standard argu- 
ments. 

Lemma 4.1. A sufficient condition for tightness of a family of probability measures supported on Sx is the 
tightness with respect to this family o/sup^gjj \K (x; lj) / F (x)\, for some continuous function F(x) which 
diverges at 00. 
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The proof is by a standard argument (which uses the Arzela-Ascoli theorem), and will be omitted here. 
Let us however note that the standard criterion for tightness of probability measures on the Skorokhod space 
requires as a second condition also the vanishing in probability of the continuity modulus wk (S) for 5 — > 0. 
That, however, is directly implied for functions in 5 a by l |21[ >. 

Lemma 4.2. For x £ R and K 6 5a , the following 'phase functional' 

$(it,Jf) := e i2 * K(x) (22) 

is continuous in K with respect to the Skorokhod topology ( restricted to S\). 

The point here is that while the evaluation function K h-> K(x) itself is not continuous, since the 
location of the jump discontinuities may change, the above phase is not affected by the location of the 
jumps when they are by integer amounts. 

4.2 The charge cocycle 

The total charge in an interval (0, u], that is: 

Q(u;u>) :=5jH - q \{j : Xj € (0,«]}|, (23) 

can be expressed as the difference: 

Q(u;u) =qK{Q;T u io)-qK{0-Lo). (24) 

In the dynamical systems terminology, Q(u\ cu) is a cocycle under the action on 5a by the group of shifts. 
That is, it transforms as: 

Q(u + s;u) = Q{u;T s uj) + Q{s-lu) . (25) 

Equation {24}, states that this cocycle is the coboundary of K, as long as the latter is well defined (a point 
which is not to be taken for granted in the infinite-volume limit). 

For background it may be of relevance to recall the following general principle. (To avoid excessive 
notation, we do not change here the notation from the specific to the general.) 

Proposition 4.3 (K. Schmidt (77)). A quantity Q(x; oS) which transforms as {25} under the action of mea- 
sure preserving transformations {T u } ue \z admits a representation similar to {24}, in terms of a measurable 
function Kq{uj), if and only if the collection of variables {Q(u; uj)} ue \z is tight. 

Tightness means in this context that the following bound holds uniformly for u £ XL 

P(|Q(u;u)| >t) < p(t) (26) 

with some p(t) which vanishes for t — ¥ oo. The less elementary part of the proposition concerns the 'only 
if direction, for which a constructive argument can be provided 1181 (as discussed also in 111). 

While the above proposition sheds light on our discussion, we shall bypass here the requirement of shift 
invariance by making use of the additional information given by Theorem |3.4| 
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Since the derivation of the enabling Theorems |3.3| and |3.4| which is presented (independently) in Sections]?] 
and [8] would take the discussion into a different arena than the one introduced above, let us first present 
their implication for our main results. 
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Proof of Theorem \3 . 1 \ - assuming Theorems \3.3\ and \3.4\ 

Existence of limits for subsequences 

Our strategy is to first discuss the question of convergence of the probability distributions at the level 
of 5a, that is of the distribution of the random particle-excess function K(x, w), which takes values in that 
Skorokhod space. Since 

sup \K(x;u>)\ < max{\K{n;t J j)\,\K{n + l;cd)|} + pW (27) 

i£[n,n+l] 

the bound of Theorem |3.3| yields: 

f( su P Iwo'm ^ - 4) - 2C ^.W) ^e-^' 1 "'^-^ 2 (28) 



rcgl 



where the upper bound is finite when t 2 A(/3, W) > 1 and vanishes, uniformly in {0, L±, L2}, for t — > 00. 
By Lemma HT1 this implies tightness of the probability measures on 5a which are induced by the Gibbs 
measures p}- p ' w ' 6 ' . Hence, for every sequence of such measures there exists a subsequence for which the 
induced probability measures on 5a converge. Since the charge configuration in any finite interval I C H 
is a continuous function of K (in the Skorokhod topology), this convergence implies also convergence of 
the corresponding point process. 

To keep the discussion simple we ignored so far (in this section) the existence of the internal degree 
of freedom yj . To incorporate that, one may consider a variant of the above argument, with a 'decorated' 
version of the function K(x; u>), for which values of the variable y are associated to the discontinuities of 
the function K. (Although this is not true for the full range of cadlag functions, for K £ 5a the collection 
of discontinuities varies continuously with K in the Skorokhod topology which is employed here.) The 
above argument applies then mutatis mutandis. 

This may be a place to note that the strategy employed here for the proof of convergence of the point 
process has its roots in the proof by A. Lenard of convergence of the Gibbs states of one-dimensional 
Coulomb systems, based on the analysis of the corresponding electric field ensemble [13 ]. For the proof of 
symmetry breaking, we employ additional arguments which were introduced in [2][T]. 

Mutual singularity of limiting measures at different values of 6; reconstruction of the phase $(w) 

The above construction of the limiting measures for the point process in T, proceeds through the con- 
struction of a limiting measure for the variable K £ 5a- In terms of this variable, the quantity we are after 
is 

$(w) = e l27TK{0 ' u ' ) . (29) 
This expression however does not suffice: for our purpose it is essential to establish that the phase $ can 
be evaluated as a measurable (and thus quasi local) function of the point configuration. That may not be 
obvious at first sight since a shift of the function by a constant: K(x) K(x) + C changes the value of $, 
without affecting the point configuration (that is the set of discontinuities of K). However, the information 
provided by Theorem |3.4| is of help here. 

Using the coboundary relation {24), for any R > 

1 r R 1 r R 

K(0;u) = -— / [K{u;w)-K(0;u)]du + - K(u;u)du 
R Jo R Jo 

1 f R 1 1 f R 

= — — / q 1 Q(u; id) du + — / K(u;uj)du 

R Jo R Jo 



£ + nu^du. (30) 
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Summing the probability bound {20\ we find that for any of the finite-volume Gibbs measures: 



(p,W,e) , , 



K(u;u)du - 1/i > S 



C 



<5 2 



(31) 



(where the insignificant l/l correction allows to relate the maximum of \K(x; within intervals of length 
A to the end-point values). It now easily follows (by considering the implications of \20\ and then applying 
the Borel Cantelli lemma) that the following limit converges almost surely with respect to any probability 
measure which is an accumulation point of the finite-volume 9 Gibbs states (the statement denoted here by 
'6 - a.s.' ) 



lim exp < i2-K 



£ I 1 

j:0<Xj<R 



R 



R_ 

■2\ 



The construction also guarantees that satisfies 



/ \ 8 — a.s. i2n 

9(uj) = e 



and by implication also the relations which were claimed in {T7J. 



(32) 



(33) 
□ 



Remarks 
1 



Extending the above considerations, one may obtain from equations $0\ and < |3 1 [ > an algorithm for 
the almost sure reconstruction, with respect to the limiting measure, of the function K from the 
location of its discontinuities. A fully deterministic quasi-local reconstruction is not possible, as the 
observation made above shows. 

The argument used above bypasses the question of translation invariance. Assuming it, or a slightly 
weaker regularity statement 1 14], the reconstruction of K (x; oj) could be done using Proposition 4.3 
through the combination of tightness of the charge cocycle, which follows from the bounds of Theo- 
rem|3.3| and the normalization of the reference level which is implied by Theorem|3.4| 



6 The one-dimensional case 

It is now left to derive the bounds which enable the above analysis. To introduce some of the ideas which 
are used in the proof in a somewhat simpler context, let us first consider the analogous question for strictly 
one-dimensional Coulomb systems. Thus we consider particles of Coulomb charge — 1 (q = 1) on a line in 
the presence of a uniform positive background charge, of density given by the Lebesgue measure (p = 1). 

We start with a bound on the probability of a uniformly large charge imbalance in such ID systems with 
appended fixed charges at the endpoints. We denote here by {x} — x — [x\ the fractional part of x 6 R. 

Lemma 6.1. For a ID jellium system of R — [r\ particles in an interval [0, r], with an 'external' charge 

7 £ N affixed at x = and charge — (7 + {r}) affixed at x — r, the Gibbs probability P satisfies, at any 
7 G (0,r): 

\ e -J(3(7 2 -l)(fl-7) 

min K(x) > 7 < , 9B ^ , , . (34) 

0<x<r V ' J e -(2fl+ 7 ln 7 ) 

Proof. We recall that in one dimension the energy of a neutral configuration is given simply by 

U(u)) = \ f K{x-Lofdx. (35) 
2 Jo 



With the affixed boundary charges, we have that K(x; lj) = 7 + x — j { i : Xi < x}\. 
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Let A be the 7?-simplex {0 < x\ < ■ ■ ■ < xr < r} of all possible particle configurations. The edges 
of A in the axes directions have length r. Consider the two subsets A + consisting of all configurations for 
which K(x; u>) > 7 uniformly in x, and A - consisting of all configurations for which: 

• There are 7 + 1 particles in unit interval [0, 1], hence K(l) = K(0) + 1 — (7 + 1) = 0). 

• There is exactly one particle per unit interval {(k — 1, k) : 2 < k < R — 7}, hence K(k) — for 
those k. 

• There are no particles in [R — 7, r], and the charge imbalance increases linearly from to 7 + {r} in 
this interval. 

For the Gibbs measure with the boundary charges as specified above, we have: 

J A+ exp[— 1/3 / K(x;oj) 2 dx] fi(duj) 



o< l x<R K ^ 'J ~ J A _ exp[- 1/3/ K(x; uj) 2 dx] fi(duj) 

< expf- inf BU(w)+ sup /3U(u) \ , (36) 

with n(duj) the Lebesgue measure on A. We shall now estimate the energy and volume factors of the RHS 
of \36) separately. 



Energy estimate: 

For uj £ A + we have (see Figure 3) 

r,, x./t + W, x e[0,l]U [R,r), 

K{x;uj) >{ D] (37) 

[7, x G [1,7?], 

thus 

2(7(cj) > A + j 2 (R- l)+B, (38) 

with A = J f) 1 (7 + {a;}) 2 dxandB = ^(7 + {x}) 2 dx. 
For uj G A - , 

f7 + {^}, x G [0, 1] U [R, r), 
\K(x)\<h, s£[l,fi-7 + l], (39) 

[7, a; £ [7?- 7 + 1,7?], 

whence 

2C/(cj) < A + (1) 2 (7? - 7) + 7 2 ( 7 - 1) + B . (40) 

Putting this together, we get the following bound on the "improvement in energy", i.e. its lowering, of 
configurations in A - in comparison to those in A + : 

inf (/-sup (7 > (7?-7)( 7 2 - 1) . (41) 

A+ A — 
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Figure 3: The function K(x) for configurations in 
A+ and A - . 

Volume factors: Applying the two-sided Stirling approximation 

\Phvn(nj e) n < n\ < V2irn(n/e) n e 1 ^ 12 

we find that 



w G A" 



and 



A.(A-) = 



< 



fl! ~ V^R(R/e) R ~ V^R 

7+1-1/12 



1 



> 



(7 + 1)! " v/27r( 7 + l)(7 + 1)'T+ 1 ' 
Combining Eqs. i ftTj l, ( [42] >, and |43j we obtain 

p («S&** (af) = 7 ) * ^y-p[-f((7 2 -D(^-7))] 

2fl+( 7 +l)ln( 7 +l) a 

xexp[-f (( 7 2 -l)(J?- 7 ))] 



g 7 + l-l/12 — *-L 2 

/3 



< e 2 «+^^xexp[-|(( 7 2 -l)(i?- 7 ))] 



(42) 
(43) 



(44) 

□ 

We now consider the ID jellium on [Li, L2], with [Li < < L2]. First, let us define two 'crossing 
events' (analogous to A of the previous lemma) which play a key role in the sequel: 

Definition 6.1. For 7, £ £ N and r > 0, let Q + {I, r) be the set of configurations satisfying 

(+a) K(0) > 3 7 and K(x) > jfor — i < x < r 

(b) K(-l) = 7 

(c) K(r—) > 7 > A'(r+). Equivalently, K(r—) = 7 + {r} arcd if( r +) = 7 + {»"} — 1. 

Definition 6.2. Le? Q = r) fee the set defined by the two conditions (b) and (c) above, but with ( +a) 
replaced with 

(a) K{x) does not equal or cross the value J for x £ [—£ + 1, R). 



THE ONE-DIMENSIONAL CASE 



Condition (a) of Q ensures that the sets Q{1, r) (I 6 N, r > 0) are disjoint. In particular, suppose H is 
the event that K(x) > 7 for at least one a>value both to the left and right of x — 0. Then [— —I + 1) and 
[i?, R + 1) are the first unit intervals to the right and left of the origin for which K(x) "crosses" 7. The 
events Q (£, r) is the decomposition of T-L into such events, thus 

H = \Jg(£,r) (45) 
Obviously, Q + C Q. We will write <j' + ' when a statement applies to both Q and <5 + . 

Theorem 6.2. 777e probability distribution of the charge imbalance, K(x), in an overall neutral ID Jellium 
in [L\, L2], with L\ < < L 2 (G Z), satisfies: 

P(|A'(0)| > 3 7 ) < c 2 e- cl ~< 3 , (46) 

for some Ci(/3), c 2 (P) > 0. 

Proof. The strategy of proof is as follows: by symmetry, it is enough to consider the case K(0) > 37 
for some fixed 7. Because of the boundary conditions K(Li) = K (L2) = (overall neutrality with no 
boundary charges), the charge imbalance must cross the line K(x) = 7 to the left and right of the origin. 
We group together configurations that have the same crossing points. We may suppose without loss of 
generality that there are no particles at integer coordinates and no two particles have the same coordinate, 
Xk £ Z and Xi < Xj for i < j. 

Denote the simplex of configurations of TV = L2 — L\ particles by 

A(N, [Li,L 2 ]) := {(xi,...,x N ) G [L 1 ,L 2 ] N : Xi < x 2 < ■ ■ ■ < x N } (47) 

and by j dto the Lebesgue integration with respect to all x-variables except the fcth one pinned at Xk = r, 
where k — \Li \ + R + 1 — 7 (recall that R — [r\ ). Let also: 

L exp(— BU (u))du) 

P(9(£, r)) ■■= r aru L ■ ( 48 > 

J A exp(— BU (u!))doJ 

Note that by the disjointness of the events Q(£, r) and ^45\ , 



I * drp(g(£ y r)) =P(«) 
e=1 Jo 



thus p(C7 (£, r))/P (H) is a probability density. 
It follows that 



P(A(0)>3 7 )< ^ / drp(a + |eMe) <supp(g + \g). (49) 



(i 



since {K(0) > 37} = U f>27 ,,>oe; + (^r). Herep(S + |£0 is using the Radon-Nikodym derivative, 
p{Q + \G) '■= p(G + ) / p{Q) where p(G + ) is defined in the fashion of {48}. 

The density p(G + | Q) is equal to the probability of (+a) occurring on [—£, r] with boundary charges as 
prescribed in (b) and (c); this is the so-called Markov property of ID Coulomb systems. Thus p(G + \ (?) is 
in the form required by Lemma[6~T) 



THE ONE-DIMENSIONAL CASE 



Note that K(0) > 37 implies £ > 27 + 1. Suppose now that 7 > y'l + ^ which in turn implies 
|/?(7 2 — 1) > 3. A little algebra gives us 

e -|/3(7 2 -i)(«+«-7) 
P (if (0) > 3 7 ) < sup , ; 



< sup e 

< e 



-c 17 2 (fl+f- 7 ) 



(50) 



for some c\{0) > 0. We can then choose c%{0) > large enough so that for all 7, P (K(0) > 37) < 



C2e 



□ 



Before moving on to the quasi- ID systems, let us give a refinement of the bound ([38j which will be 
useful for proving an analogous result to Theorem |6.2| for the strip. 

Lemma 6.3 (Improved energy bound). Let u G Q + . Ifn k (w) is the number of particles in [k, k + 1] for 

k — —I, R — 1 and 71r(w) is the number of particles in [R, r), then n—t(w) = and 



K(x;u J ) 2 dx>" f 2 (£ + r)+ ^ (jn k (uj) 2 + n k (cj) 3 /S). 



(51) 



Proof, n-i = follows from K(—£+l) = 7+ 1 — n-i > 7 and the integrality of the involved quantities. 
More generally, n k < k + £ and J^^L-; n k < m + £ for all m < R. The energy estimate follows from 
the observation that any non-zero n k adds at least a triangle (see Figure 4) with side length n k to the field 

K(x),x < k, on top of 7: 

K{k) = 7+ (k + £) - > -y + n k . 

3<k 



This gives additional factors in the energy 



2'yxdx + / x dx=^fnk{oj) + n k (uj) /3. 
11 Jo 
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(52) 
□ 



Figure 4: Excesses in the energy of ui £ Q + , as estimated in \5i\ . 
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7 Tightness bounds 

We now turn to the proof of Theorem |3.3| for jellium in the tube T. The strategy is as demonstrated 
in the ID case. However, the energy estimates are complicated by the V2-interaction. We denote by 
V 2 (oji , 0J2) the sum of the V2 -interaction terms of pairs of distinct particles, one in wi and the other in 012, 
and by V 2 (u) := V 2 (w, w) such a sum for a single configuration, omitting self interactions. Following the 
decomposition (|7}, the energy of a configuration uj = [z\, zjv) is a sum of a ID part and a short-range 
interaction energy. 

U(u)= ^ \K(x;u)\ 2 dx + q 2 Y, V 2 (y j ,y k ;\x j -x k \) = U 1 (u) + V 2 (u). (53) 

l<j<fe<JV 

Given £ > 0, r > 0, the tube T splits into subsystems 

Y = Y t , r ■■= [-t, r] x D, CTZ := T\Y. (54) 
and the total energy naturally decomposes into 

U(w) = Ucn(^cn) + Uy(uy) + V 2 (lo y ,ujctz) ■ (55) 

Note that in one dimension V2 = so that the last term on the RHS of ^55\ vanishes. This yields the 
Markov property for the function K(x; u) in one-dimensional Coulomb systems (where it plays the role of 
the electric field). 

In the current setting, we must control 

Vifwv) + V2(wy,W£7j) ■ (56) 
Following are two useful perspectives on the arguments which are presented below. 

An alternative viewpoint: Replacing particles with rods To analyze V2 -interactions we can 
employ the procedure of replacing a particle at with a rod, or in dimension d > 3 a D-shaped charge 
slab, of vanishing thickness = {z : x — x k }, in which the charge (— q) is uniformly distributed with 
charge density —qW~ 1 dy. By {3J the 'rod-particle' or 'rod-rod' interaction is completely free of the term 
V2, and thus when the particles in Y are replaced by rods, the rod's distribution conditioned on the rest of 
the system is identical to that of a strictly one-dimensional system. In this setting, one will broaden the 
notion of a configuration of particles to a generalized configuration of particles and rods. The x-values will 
still be ordered, and the replacement z k i-» corresponds to the map 

u = (zi, . . . , zjv) h-> u; {xk} = (zi,. . . , z fc _i,x fe , z k ,. . . , z N ) . (57) 

Each replacement changes the energy by an amount equal to (56}: Thus, denoting by uj 1— > u/" (generalized) 
configuration in which all particles in Y were replaced by rods, we have 

U(oj) - U(lj y ) = V 2 {cjy) + V 2 {m,ucTi)- (58) 

Therefore there are two equivalent points of view for the estimates of the next section: we can either 
think of V 2 -energy estimates for a standard system of point particles, or we can think of the estimate as 
bounding the replacement effect on the energy incurred when all particles in Y are replaced by rods. 
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.Tellium as a system of unbounded spins One may also regard the system discussed here as one of 
unbounded spins, with additional internal degrees of freedom. There is however, a significant difference, 
which does not make the existing bounds for such systems with 'superstable interactions' 1151 [T6l [T2l 
directly applicable. The "spins" are the charge imbalances at integer x coordinates, K(k\ u>), k G Z. The 
internal degrees of freedom are the exact locations, including y coordinates, of particles in [k, k + 1) X D. 

The V2 -interaction between two disjoint regions X and Y is bounded below by -q 2 g(dmt(X, Y))n(X)n{Y), 
with n(X) the number of particles in X. This is reminiscent of lower regularity as defined in 1151 . Upper 
regularity, i.e., a (pointwise) converse inequality, does not hold since V2 = 00 is possible. However, a 
suitable substitute can be obtained through the zero mean of V2 over y (applying the Jensen inequality). 

The representation \52>) together with V2{cu) > —q 2 X]j<fc nk(u))nj(ui)g{k — j — 1) suggests that an 
inequality of the form 

U(ui) > ^[Afsrffc; cj) 2 - B] (59) 

k 

for suitable A, B > might hold, i.e., the system of spins might be superstable. 

This is really close to situations considered in 115111 61 - In fact, the proof of Theorem [33] is close in 
spirit to the proofs in those papers. One complication arising for jellium is that we cannot simply decrease 
one spin K(k; oj) without affecting other spins: point charges cannot be removed without affecting an 
infinite change in the energy, but they can be moved from one place to another. In other words, the overall 
neutrality fixes the total number of particles, in contrast to the grand canonical setting of 1151 . The procedure 
employed here circumvents this difficulty: as in the ID case, we look not at an individual large spin but at a 
selected subsystem [—£, r] x D = Y, inside of which charges are rearranged leaving the charge distribution 
outside unchanged. A technical difficulty here, shared with the systems of I15II16I . is that making some 
spin smaller may actually increase its interaction energy with other spins. 

In the rest of the section, we shall write P instead of Mjf^'^p an d without loss of generality, we will 
assume 8 = 0. 

7.1 VS>-eiiergy estimates 

We would like bounds on what a replacement does to the total energy and to the Boltzmann weight of a 
configuration oj. The first step towards this end is the zero average property {5). Together with Jensen's 
inequality applied to the exponential function, it yields 

e - Ww fl*>) < 1 r e -m») dyh . (60) 

" Jo 

Our goal for the rest of this subsection will be to prove a sort of converse to l |60| >. 
Let I € N, r > and Y C T as in {54). Let R = [r/X\ . We divide Y into cells 

Y k = [kX, (fc + l)A) X D for k = -£,..., R- 1 and Y R = [RX, r] x D, (61) 

and denote by n^iuj) the number of particles in 

We start by noting that V2 is "lower regular" in the sense that it can be bounded below in terms of the 
(w-dependent) particle numbers nk(w). In this statement a role is played by: 

Lemma 7.1 (Lower bound on V^(wy)). For all u>, the Vi-energy satisfies 

R 

V2 (wy) > — qC n k (uJY) 2 (62) 

k = -l 



for some constant C = C(W) > 0. 
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Proof. By the assumption which was made in {9}, the interaction V2 between particles, in cells Y kl , Yfc 2 
(not necessarily distinct) which are distance DA > from each other, is bounded below by 

- q 2 g{D\) n kl (uj) n k2 (ui) , (63) 

with g(u) the non-negative monotone decreasing function which appears there. The integrability assump- 
tion on g implies that the sum over integer multiples of A is finite, X^dgn ffC^^) < 00 • 

Let ko be the value of k such that nfc(w) < n k() (uj) for all k = —£, . . . , R. The interaction between 
particles in Y k (or Yr = [i?A, r] x D if k = R) with themselves and with others is lower bounded by 

fc -i R 
-q 2 g(0)n ko (u) 2 /2-q 2 ^ g([k -k-l]X)n ko (uj)n k (ui)-q 2 ^ g([k-k - 1] \)n k(> (u))n k (bj) 

k=-i fc=fc + l 

> -I 2 ff(0) + 2 f>(DA) j n k0 (u) 2 =: ~q 2 Cn ko (uj) 2 . (64) 

Next, we can find fci maximizing n k (uj) among k 7^ ko, and estimate interactions between the remaining 
cells, k — —£, . . . , R, k 7^ ko, in a strictly analogous way. Iterating the procedure until all cells have been 
taken care of, yields the lemma. □ 

Next, we need to control the interaction between Y and the rest of the tube (CTZ). Obviously, as in {63), 
the interaction between particles in, e.g., Y-e and the left subsystem may be bounded below by 

Li -e-i 

- q 2 22 g(DX) n-i («) n_£_D_i(w) . (65) 

D=0 

However, for long tubes, with L\ » 1, this bound can in principle become very negative. That happens 
in case there is an accumulation of particles of C1Z close to Y. We will need to control the frequency of 
occurrence of such "irregular" configurations. 

Let 76 N. Let 5 C Jibe defined as in Section|6]with suitable y degrees of freedom added. Recall that 

for oj g Q, 

K{-£; u) = 7, K(r-;u>) > 7 > K{r+; w). (66) 

We will call a configuration cj G G regular if its charge imbalance is well-behaved outside Y in the sense 
that 

\K(x;oj)\ <7 + dist(a:, [-t,R])forx$ [-£,r]. (67) 
The set of regular configurations will be denoted here Q Icg C Q. 

Remark: On Q, we automatically have that K(x; lj) > —7 — | x + l\ for x < —£ and that K(x; tj) < 
7 + I a; — {r}\ for x > R (see the regions corresponding to the dotted triangular regions in Figure 5), 
thus the condition \67\ really amounts to restricting K(x; u) from entering the shaded triangular regions of 
Figure 5 below. 
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Figure 5: Regular and irregular configurations 
[j G 6, in the sense of {67\. For regular 
I configurations the graph of K does not enter 
the shaded region. 

Within [-£, R + 1], the function if can 
cross 7 only in the first and last unit inter- 
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For configurations in C/reg we can lower bound the V2 -interaction and the total energy change due to a 
replacement: 

Lemma 7.2 (Lower bound for VifW, wen))- For all ui £ <5 rC g, ffe Wenergy satisfies 

Viiwv.tJcn) > —<lC"i max ndw) (68) 

-l<k<R 

for some constant C = C(W) > 0. 

Proof. For fc = —t, R let rifc(w) be the particle numbers introduced above. We extend the definition to 
k £ Z, Li < k\ < 1/2 — A in the obvious way. Call n r (uj) the number of particles in (r, (R + 1)A) x D. 
Using the monotonicity of g(u), the interaction between a particle z S Y with the particles from w in 
s < — £ is lower bounded by 



q 2 j([Hi + «]A)n_,_ H ( W ) > -2q 2 Y,9([k + x + £]X). 



(69) 



The inequality is obtained as follows: sum by parts in order to rewrite the sum with differences of g and 
sums of particle numbers, use that g is decreasing and oj regular, and sum by parts again. Thus the total 
interaction of all particles in Y with the left substrip is bounded by 



-2q 2 ^ n p (uj)^2g([p + £ + k]\) > -q 2 { max R n k (cu) j ^ 2(k + l)g(kX) 



The interaction with the right tube x > r can be bounded in a similar way. An additional summand 7 will 
appear because for regular configurations, there may be as many as 27 particles accumulated near x — r 
(see Fig. 4 above). □ 



Lemma 7.3 (Coupling lemma). For all ujy (a configuration of particles inside Y), 



L 



-PV(w) 



Ch >/2 +nM BI (») 3/! +7»„,x(«)) 



(70) 



for some constant C = C(f3,W) > and n max (^) = max_KKij nfc(w). The integrals are over 
configurations wen where 

W '.= Wen U CJy G S(reg)- 
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In somewhat loose notation, Eq. l |70| l may be thought of as a lower bound for the probability of regular 
configurations, given that they are in Q and that the configuration inside Y is ojy 

P(<?r C g | Q,m) > exp(-C(7 3/2 + n ma x(^Y) 3/2 + 7n max (a;Y))) (71) 
In this sense Lemma [73] says that regular configurations have high enough probability. We defer the 



Using Lemmas 7.1 7.2 and 7.3 we can now formulate the "converse" to Eq. l|60|l. Note that U(ui ) is 



the energy of the system with all particles in Y replaced with rods. 
Lemma 7.4 (Replacement lemma). For all loy, 

-W^dum < e c ^ 3/2 +~^>+£?=-* / e-W^ducn, (72) 

e Jq 

for some constant C = C(j3, W). 

Thus, on average, the Boltzmann weight before replacement is smaller than the Boltzmann weight after 
replacement, up to some (controllable) function of wy. 

7.2 Proof of Theorem @T 

The strategy of proof is exactly the same as in the ID case. Let 7, — L\ , L2 G AN. We want to estimate the 
probability of the event K(Q; uj) > 37. Let (£, r) be events defined as in Section]^ Note that since the 
events were defined in terms of charge imbalances which depend on ^-coordinates solely, we can extend the 
definitions to the quasi- ID case by adding the y degrees of freedom. The definition of conditional densities 
p((5^ + '), p(Q^\G) is extended in the natural way as well. By the same argument as in Section]^ 

P(ff(0;w) > 3 7 ) < supp(g + (£,r)\g(£,r)) (73) 

On CJrcg, Lemmas 1 7 . 1 1 and |T2"| give us control over the V2 -interactions, enough so that the "ID-portion" 
of the energy implies the desired exponential decay of the RHS of l |73[ ). However, to deduce the exponential 
decay for all of Q, we need Lemma |7"!4] (which combines Lemmas |7.1|7.2| and |7.3} . 

Recall that U(lo y ) = Uy(ljy) + Uck(uctz) is the energy of the system with particles in Y replaced 
with rods. By Lemma|7~4] 



e -pV{u) < / e -/Jir(J) e C( 7 3/2 +7«„ax("Y)+Ef = - ( »t ("y) 2 ) dCj (74) 

where we recall that dw means Lebesgue measure of all variables except the a>variable pinned at x = r, 
and by Jensen and Eq. ([S} 

/ exp(-pU(u)) du) > f exp(-/3C/(tj Y )) dCj. 

It follows that 



p{Q + \9) < 



f g+ e X p(-/3[/(u; ¥ )) e X p[C( 7 3/2 + jn^oj) + J2k=-e n k {ujY) 2 ]duj 



Note that this upper bound is independent of L\ and L2. 

Since in U (^ ¥ ) there are no interactions between Y and CTZ, we have again a Markov property, and the 
estimates reduce to estimates for the finite tube Y. We can proceed as in Lemma [6T| using Lemma |63| in 
order to take care of the extra factor exp[C"(7 3,/2 + 7n max (cjY) + Ysk=-e n k(uY) ]• 
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7.3 Regular configurations have positive probability 

In this section we prove Lemma [73] The strategy of proof is to define a map 

G — > Greg, w i— > UJ (75) 

by shifting some of the particles in C and 1Z (the left and right substrips in the complement of Y) away 
from Y, and to substitute jg duj by J g do/. Complications arise because the map is not one-to-one and 
has a non-trivial Jacobian, resulting in a compression of phase space, i.e., entropy loss. However regular 
configurations have lower (ID) energy, and this energetic improvement is enough to compensate for the 
entropy loss. 



Remark. The astute reader will notice the parallel between the above map and the comparison (in the ID 
case) of A + and A~. Configurations in A - are energetically more favorable but have smaller Lebesgue 
volume (i(A~). 

Fix 7, £ G N and r > 0. Let R — [r/X\ . For uj G G, denote ujy , ujc and uj-r. the projections onto 
Y, £ and 1Z. To simplify matters, let us focus on 1Z and pretend first that I = L\ (which really cannot 
happen since K(—£; uj) = 7 > = K(L 1 ;uj)). There are iV-^ = L2 ~ R — 1 + 7 particles in 1Z. Write 
ujtz = {zi, ZL 2 -R-i+t} with particles labeled from left to right, Xj < Xj+i. We observe that if for all 

k 

x k - r > (k - 7 - l)A/2, (76) 

then uj € Greg- We will call particles regular if they satisfy Eq. j76[ >, irregular otherwise. Thus regularity 
of all particles implies regularity of the configuration. The first 7 + 1 particles are always regular. 

For uj £ G, let uj 1 G G be such that uj' y = ujy (no changes in Y) and uj^ the collection of points {z' k } 
with 

, , (xk, regular particle, 

[ (R + k - 7) A + (fc _^_ 1)/2 , irregular. 

The idea behind the map is the following: for fixed charge imbalance K(r) ~ 7, the ID part of the energy 
in 1Z is minimal if 7 particles accumulate near x = r and the remaining L2 — R — 1 particles occupy each 
one of the "equilibrium" cells Y r+i, ...,Yl 2 -i- The map uj i-> uj' simply shifts an irregular particle closer 
to its equilibrium position, thereby decreasing the (ID) energy. 

Lemma 7.5 (Energy estimates). Let Ii TT (uj) C {7 + 2, ...,N-jz} be the set of irregular particle labels. The 
map G 3 UJ <— > uj' G Srcg decreases the ID-energy by an amount which is at least 

U\u)-U\u') > ^ £ (*-7) a -l , (78) 



(«0 



Furthermore, 



V 2 (ujy,ujiz) > -g 2 Cn max H(|/ irr (o;)| + 7+ 1) (79) 
V 2 (uJn)>-q 2 C + ( 80 > 

k£I irT (ui) 

for some constant C — C(W) > 0. 

Note that we do not evaluate Vi(w) — V2(uj'). Instead we estimate directly V2(uj). Jensen's inequality 
will take care of V2(u/). 



TIGHTNESS BOUNDS 



Proof of Lemma [Z5| The ID-energy term: We shift irregular particles successively by the order of their 
subscript, starting with the highest one. In this way one obtains a sequence ufi 1 ' starting at a/ ' = ui. The 
shift x k — > x' k increases the particle imbalance in [xk,x' k ) by 1 and leaves it unchanged elsewhere. This 
decreases the energy by an amount 

„2 rx'u 



u 1 ^) - c/ 1 ^" -1 ') = 



2W 

w 



\K(x;J n) f - (K{x-J n) ) + l) }dx 



(81) 



where we have used that the integrand is bounded above by an affine function with slope 1/A, end value 
< 1/2, integrated over an interval of length x' k — x k > (k — j)X/2. The sum of these bounds yields JTSJ. 
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Figure 6: A possible energy improvement for an irreg- 
ular [j £ 8, which is obtained through the displace- 
ment of a particle from x = a to x = b. The modified 
K(x) is described by the dotted line. 

V2-energy inside 1Z: The individual terms Vi are bounded below by — q 2 g(0). Since there are 
Xlif=i ' 1 & pairs of irregular particles, we get 



XI V/ 2(yj>yk;\xj - x k \) > -q 2 g(0) ^2 k. 



(82) 



j<k 



The interaction between a given irregular particle z k and the regular particles in 1Z is lower bounded by 



- q 2 ^2g(\ Xj - Xk \\)>-(k-l)q 2 g(0)-q 2 £ g((j - k)X/2). 

3 reg j>fc+l 

It follows that 

oo 

£ V a (y s , y k ; \x 3 - x k \) > -g 2 X(s(°) fc + J2 9(n\/2)) 

3 rc g k irr. ri — 1 

k irr 

Vi -interaction between Y and 1Z: Eq. \19\ holds provided C is large enough so that 

oo oo 

^^ 9 (mA+[n/2]A)<C. 

m— n— 

Lemma 7.6. For all a;, 7 and £, r, Li, L2, 



(83) 



(84) 



(85) 
□ 



/ exp(-@U((jj))dyi ■ ■ ■ dy Nlk 
Jo N n 



< C" exp(C"( 7 3/2 + n max H 3/2 + 7nm„(^))) exp -C ^ (fc - 7 ) 2 ) 



/ exp(-/3U(u'))dyi---dy NK , 
Jn N n 



(86) 



TIGHTNESS BOUNDS 



for some constants C = C(/3, W), C = C'(/3, W), C" = C"(f3, W). 
Proof. 



J 



-PU{u) dy < e - / 9[!7 1 (^)-(7 1 ( lJ ')+mf I) (V 2 (^ TC ) + V 2 (^ Y ,^ TC ))] / eT^^dy 



< e -nu 1 (^)-u 1 {u 1 ')+inf v (v 2 (u 1 - R )+v 2 (u lY ,u 1 - R ))] / e -P u{ul ">dy (87) 



We might think of this inequality as a result of a three step procedure: (a) Replace particles with rods. The 
replacement error is quantified by Eqs. (80) and (79) . (b) Shift irregular rods away from Y. The energy 
decrease is given by Eq. l |78[ l. (c) Replace rods with particle averages using Jensen's inequality. 
Next, by Lemma [73| letting 77 = k — 7 



U^lu) - U\J) + in{{V2{LUn) + Vi(uy, CJ-fl)) 
y 

2 \ 

> [|^7^ 2 - C(7 + + "max + 1)J -Cn max (cj)(7+1) (88) 

V- (l + 7) S-Tirr 

The Crj is dominated by rf except for small r\. The 7 term inside the sum is dominated by r/ 2 except for 
rj 2 < const 7, whence a negative contribution of the order of r y 3 ^ 2 , and similarly for n max (u;). Thus the 
RHS of ([88} is greater than 

Ci ^-7) 2 -C , 2 (7 3/2 +n max H 3/2 +7n max H) 
fcei lrr (w) 

for suitable constants Ci = Ci(/3, W), C2 = Ci(fi, W) (at this point the constants depend on /3 only 
through their dependence on q; see the remarks concerning dimensionless parameters in Section |2~2} . □ 

We are now equipped for the proof of LemmaPO] 



Proof of Lemma [ZJ| Let 7, £, r be fixed. As in previous proofs, we will deal only with the region TZ (C is 
dealt with analogously). 

Define the map ui 1— > <J from Q to Q rcg as above. Let I C {7 + 2, N-jz}. Conditioned on the event 
that Jiir(w) = /, the regularizing map is injective with the Jacobian 



duj' 



k 



n fc^-r - ( - e m(* - ^ - 1)/2]) 

fee/ ' \ fee/ / 



(89) 



Thus, for suitable constants, 



Ji iTX (u)=i Jg r , 

with 



e-^ ( -''d^, (90) 



A(7) := /3C^> - 7 ) 2 - PC'(y 3/2 + n^ x + 7 n max ) - C ^ ln[(fc - 7 - l)/2] , (91) 

fee/ fee/ 

where the last sum comes from the Jacobian. We obtain a lower bound on A(I) by dropping the last sum 
and adjusting the constants C and C' . It follows that 

00 

Y exp(-A(/)) < e c 'h 3/2+ "- +1 "»») Y[(l + e- Cfe2 ). (92) 

/C{ 7 +2,...,iV TC } fe=l 
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Summing Eq. {90) over / yields Lemma |7"3j up to the simplification I = L\. 

In the general case I < Li, the regularizing map is extended by defining its action on particles in C in 
a similar way: particles too close to Y get shifted away. The energy estimates of Lemma [73] extend in the 
natural way. From there one can proceed as in the fictitious case I = L\. □ 



8 Convergence of the volume-averages of K(x; uj) 

In this section, we prove Theorem |3.4| In order to acclimate ourselves, let us first prove the ID case via 
a shift coupling for fixed 5 > 0, — Li,I/2 G N, and r as in Theorem |3.4| As before, without loss of 
generality, we assume 9 — 0. 

Recall that Xi denotes the position of the zth particle with i < j implying that Xi < Xj almost surely. 
The coupling a; M> <!; is defined by moving the particles in the following manner: 

• If Xi < 28 then Xi = Xi/2 

• If 2S < Xi < r then Xi = Xi — S 

By symmetry, we may think of Theorem |3.4| as a statement concerning the exponential decay (in r) of 
the probability of 

rr 

A s := {uj : / K(x-w) dx > 8r}. 



To prove this exponential decay, we use the observation that using the above coupling, for 25 < x < r, 

K(x;u)-5 = K{x-8;u). (93) 

Proposition 8.1. Consider a one-dimensional Coulomb system. For all 5 > and c < 1, there is a 
constant ro(5, c) such that r > ro implies 

V(A S ) < e- cfis2r . (94) 

Proof. If we define the set 

T = {uj : \K(k; w)| < r 2 ' 5 - 1 for each integer k £ [0, r]}, (95) 
then by Theorem|6.2|we have 



P(^ c ) < (r + l)c2(P)e- clWr . 

Note that if u G T then \K(x;cj)\ < r 2/5 and \K(x;Q)\ < r 2/5 + 5 for all x G [0,r], Thus for 
uj G A s n 7: 

K{x\ u) dx>5r- 2Sr 2/5 , (96) 



J2S 



and using |93} 



K(x;uj) 2 -K(x;u)fdx > / (2SK(x; uj) - 5 2 ) dx - / K(x; Co) 2 dx. 

J2S J [0,6]U[r-6,r] 

> 25(Sr - 2Sr 2/5 ) -S 2 (r- 2S) - [ K(x; Cj) 2 dx 

J [0,S]U[r-5,r] 

>8 2 {r + 26- 4r 2/5 ) - 2S(r 4/5 + 25r 2/ '° + 5 2 ) 

= 5 2 (r - 8r 2/5 - 25- 1 r i/5 ) (97) 

To get a bound on the probability we now only need a bound on the change in the Jacobian or "volume- 
factor" caused by going from uj to uj. The only place where volume is lost is in the first step of the coupling 
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when [0, 25] is mapped to [0, 5]. For u £ T there are at most 2r 2 ^ 5 + 25 particles in [0, 25], thus the 
volume factor is bounded above by e' 2r +2S ^ log 2 . Altogether we have 

P(A 5 n^ c )+P(^n^)<(r+i) C2e - cl '' 6/ %^^ 

(98) 

from which {94) readily follows. □ 

For the general quasi-lD setting, the coupling uj h-> uj maps a;i t-> Xi just as above and leaves the 
2/j's unchanged (j/j = We need to adjust the proof of Proposition [8T| in two places. The first trivial 
adjustment is to bound the probability of T c with Theorem 1 3 . 3 1 (replacing the ID bound of Theorem |6.2[ l. 
The non-trivial adjustment is to show that on T, any increase in the short-range V 2 -interaction energy 
caused by the coupling is insignificant compared to the decrease in ID-energy given in l|97|>. 



7.1 



Proof of Theor em\3.4\ Let T be as in (|95J and define the 'regular set' J- rBg similarly to the definition of 
with r 2 / 5 playing the role of 7. By Lemma ' 



7.3 



we have that 



P(J- rcg ) >e~ cr4/5 ¥(T) (99) 

for some c > 0. We now condition on J- TCg and bound J exp(— dy j J exp(— pViiCj)) dy on 
this event. 

Recall (from Section|7]l that Y 0i i, = [a, b] x D and that uj Ya - b is ui with all particles in Y a ,j, replaced 
with rods. If there are k particles in the region Yo,s (for Co) and the region ¥0,25 (for lo), and y denotes the 
vector of ^-coordinates of those particles, then for every u G J 

/ exp(-/3V 2 (w))d2/ 

<exp(-^ inf [V 2 {lu)-V2(Cj Yo - s )]) [ e X p(-/3V 2 (Co Yo - s ))dy 

< expf-/J inf \Va(u) - F 2 (w Y °- 2S )] - P[V 2 (uj Yo ^) - V 2 (Cj y °^)]) f exp(-^y 2 (a;))dy, 

(100) 

where we have compared the integrals in the last two lines using Jensen's inequality applied to rod replace- 
ments. We have slightly abused notation in that the infimum over y £ D fc is at different rr-values for the 
coupled configurations Co and uj. 

Note that the coupling does not affect V 2 -interactions between two particles that are both in ¥2,5,,- (pre- 
coupling). Therefore, on the event J- TCS , the affected V 2 -interaction can be lower bounded by 

V 2 (uj) - V 2 (UJ C TZ) - V 2 (lV Y \Y 2S ) = V 2 (u}Y,U}£Tl) + V 2 (uJY 0t5 ) + V 2 (LJY . s ,U}Y\Y _J > -Cr 4/5 . (101) 

tied on J- Teg , 

V 2 (io Y °- 2S )-V 2 (uj)<Cr i/5 (102) 



If we can now show that conditioned on J- reg , 



and 

V 2 {Cj y °- s ) - V 2 (uj y °- 2S ) < Cr 4/5 , (103) 

then the ID decrease in energy caused by the coupling (see |97j) overwhelms the coefficients of the right- 
hand sides of both < | 1 00| > and l |99| l which would prove the theorem. Let us now show l |102} and ( | 1 03 1 > - 

To bound V 2 {uj V() - 2S ) - V 2 (u) on .F rog , note that there are at most 2r 2/5 + 25 particles in Y , : 
Replacing the role of max rife with 2r 2//5 + 5 and the role of 7 with r 2 ^ 5 in Lemma 



7.2 



we obtain 



V 2 (oj y °- 25 ) - V 2 (uj) < C{r 4/5 + 5r 2/5 ) (104) 
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for some C > 0. 

To bound V 2 (w Y ° • a ) — V 2 (<jJ Y ° ■ 2 5 ) , first recall that ijy r and ujy^ denote the sets of particles of uj that 
are in Yo, r and its complement, respectively. Using l |58[ >, since the V2 -interactions within Yo, r are the same 
forbothd; Yo ' 5 and uj Yo - 2S , we have 



By (TO) we have that for each S > there exists C(<5) > such that 

^2(^1, 22) < Cg(\xi — x 2 \) 
whenever \x\ — x 2 \ > 8. The condition \xi — x 2 \ > S is satisfied by 21 £ oDy° '* and z 2 £ lip s 



(105) 



following the proof of Lemma 7.2 we get that V 2 (u> 



so again 

V 2 (uj) < Cr 4/5 for some C > 0. □ 



9 Discussion 

The results presented here confirm a conjecture which was stated in (TJ. Let us point out some related 
questions which were not addressed in this work. 

1. It may be worth stressing that the length A = (pW)" 1 , with which translation symmetry breaking 
is proven here for shifts by u ^ XL, does not correspond to the interparticle spacing. The difference 
between the two length scales shows up only for values greater than 1 of the dimensionless parameter 
M^/A^ -1 ', but it becomes very pronounced when W ^> \^ d ~ 1 ' . In particular, this means that the 
symmetry breaking proven here does not correspond to the phenomenon of Wigner lattice. 

2. A question which was not addressed is whether the symmetry breaking stops at the level which is 
proven here, e.g. whether for each 9 the Gibbs measures pi^'^f 1 have a unique limit, for n\,n 2 — > 
00. The answer would be negative if, for instance, for at least certain values of W the system is in a 
lattice like state wrapped on the cylinder. 

3. Let us express here the conjecture that the answer to the above question is positive. More explicitly: 
we expect that for each 6 there is a unique limiting state. Furthermore we expect this state to be not 
only invariant under the shift T\ (which would be implied by the uniqueness) but also ergodic with 
respect to it. In particular this would mean that the states do not admit any further cyclic decomposi- 
tion; i.e., the symmetry breaking stops at the level which is proven here. Can that be shown in a brief 
argument? 

4. Adding a comment to the above question: in case of the cylinder we do not expect the symmetries of 
rotation in the compactified dimensions to be broken. 

5. As is the case for point processes which are not too singular, the ^-states discussed can be charac- 
terized through their correlation functions, {p n (z 1 , . . . , z n )}. Thus, the state's non-invariance under 
shifts implies that at least some of these correlation functions (if not all) are not shift invariant. In (9) 
it was shown that in narrow enough strips translation symmetry breaking occurs already at the level 
of the one-point density function. Does such a statement extend to the entire regime covered by the 
non-perturbative argument presented here? 

6. Other examples of particles with Coulomb potential include two component systems, of particles of 
charges ±g (and not necessarily equal masses). For such systems on a line, there is no translation 
symmetry breaking, but there is a related phenomenon of phase non-uniqueness j2j- In essence, 
in one dimension the fractional part of a charge placed at the boundary cannot be screened by any 
readjustment of the integer charges. Does this phenomenon also persist to the quasi one-dimensional 
Coulomb systems? We conjecture that the answer is affirmative even though the rigidity of the ID 
Coulomb interaction is somewhat 'softened' in the quasi one-dimensional extension. 
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